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Abstract 

We investigate nuclear matter properties in the relativistic Brueckner approach. 
The in-medium on-shell T-matrix is represented covariantly by five Lorentz invari- 
ant amplitudes from which we deduce directly the nucleon self-energy. We discuss 
the ambiguities of this approach and the failure of previously used covariant rep- 
resentations in reproducing the nucleon self-energies on the Hartree-Fock level. To 
enforce correct Hartree-Fock results we develop a subtraction scheme which treats 
the bare nucleon-nucleon potential exactly in accordance to the different types of 
meson exchanges. For the remaining ladder kernel, which contains the higher or- 
der correlations, we employ then two different covariant representations in order to 
study the uncertainty inherent in the approach. The nuclear matter bulk properties 
are only slightly sensitive on the explicit representation used for the kernel. How- 
ever, we obtain new Coester lines for the various Bonn potentials which are shifted 
towards the empirical region of saturation. In addition the nuclear equation-of-state 
turns out to be significantly softer in the new approach. 
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1 Introduction 



The investigation of nuclear matter properties within the relativistic Dirac- 
Brueckner-Hartree-Fock (DBHF) approach [1-6] remains a fundamental topic 
in theoretical nuclear structure studies. Compared to non-relativistic approaches 
the relativistic DBHF treatment turned out to be a major step forward in 
the explanation of the saturation mechanism of nuclear matter. The satura- 
tion points obtained for non-relativistic calculations, throughout all possible 
choices of different nucleon-nucleon interactions, are located on the so called 
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'Coester line' [7] which does not meet the empirical saturation region. Us- 
ing modern nucleon-nucleon interactions of the one-boson exchange type [8,9] 
the relativistic calculations also reveal such Coester lines which are, however, 
significantly shifted towards the empirical region [3]. 

On the other hand, many details of the relativistic theory are still not fully 
resolved. In particular, the precise form of the nucleon self-energy, i.e. the 
magnitude and the momentum dependence of the scalar and vector self-energy 
components are a question of current debate [3-5,10,11]. Since the self-energy 
describes the dressing of the particles inside the medium and thus determines 
the relativistic mean field this fact states a severe problem. Different techniques 
to handle the DBHF problem can lead to significantly different results [3-6,11]. 
In a recent work [5] we found that the momentum dependence of the nucleon 
self-energy is dominated by the one-pion exchange contribution which accounts 
for the nuclear tensor force. 

Unfortunately, the treatment of the 7rNN vertex and the corresponding self- 
energy contributions is closely connected to a severe ambiguity in the T-matrix 
representation [5,10]. The DBHF approach starts from a realistic nucleon- 
nucleon potential of the one-boson exchange type, i.e. the Bonn potentials [12]. 
As for the free two-body scattering problem, anti-particle states are neglected, 
and thus one works exclusively with positive energy states. Hence, a direct 
determination of the nucleon self-energy operator is not possible since not all 
matrix elements of this operator are known. Horowitz and Serot have therefore 
developed a projection technique to determine the scalar and vector self-energy 
components from the in-medium T-matrix [1]. In this approach the T-matrix is 
represented covariantly by Dirac operators and Lorentz invariant amplitudes 
where the latter are determined from the positive-energy on-shell T-matrix 
elements. 

The whole problem arises from the fact mentioned above, namely that one 
does not include negative energy states and therefore neglects the excitation 
of anti-nucleons. The inclusion of negative energy excitations with 4 states for 
each spinor yields 4 4 = 256 types of two-body matrix elements concerning their 
spinor structure. Symmetry arguments reduce this to 44 for on-shell particles. 
[13]. If one takes now only positive energy solutions into account this reduces 
to 2 4 = 16 two-body matrix elements. Considering in addition only on-shell 
matrix elements the number of independent matrix elements can be further 
reduced by symmetry arguments down to 5. Thus, all on-shell two-body matrix 
elements can be expanded into five Lorentz invariants. But these five invariants 
are not unique since the Dirac matrices involve always also negative energy 
states and thus a decomposition of the one-body nucleon-nucleon potential 
into a Lorentz scalar and a Lorentz vector contribution depends on the choice 
of these five Lorentz invariants mentioned above. The best choice would be 
to separate completely the negative energy Dirac states. But since this is not 
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possible, there is not a unique but only an 'optimal choice'. The topic of this 
paper is the form of this 'optimal choice' of the five invariants. 

Thus, as discussed in [5], various covariant representations of the T- matrix 
exist which all reproduce identically the on-shell T-matrix elements but lead 
to rather different results for the nucleon self-energy in the present formalism. 
This ambiguity was found to arise substantially from the treatment of the 
pion exchange part, in particular by the way how to take the pseudo- vector 
nature of the 7rNN vertex into account. In realistic nucleon- nucleon potentials 
the pion is usually described by a pseudo- vector coupling [12]. A pseudo- 
vector 7rNN vertex is also predicted by the non-linear a-model based on chiral 
symmetry [14]. Furthermore, a pseudo- vector 7rNN vertex has the advantage 
that it is consistent with neglection of negative energy states while a pseudo- 
scalar coupling connects very strongly positive and negative energy states. 

Following this argumentation in several works [2,4,10,15] the so called 'pseudo- 
vector choice' to the T-matrix representation was applied. This means to sim- 
ply replace the pseudo-scalar by a pseudo-vector covariant while keeping the 
corresponding amplitudes unchanged. However, this procedure is not well de- 
fined since equivalent ps representations exist which lead to non-equivalent pv 
representations. A strong momentum dependence of the nucleon self-energy 
emerges, e.g., in the 'conventional' pv representation as it was applied by Sehn 
et al. [4]. In Ref. [5] we addressed the problem with the above ambiguity, in 
particular with respect to the determination of the nuclear self-energy compo- 
nents. There we also proposed a 'complete' pv representation of the T-matrix 
which results in a much weaker momentum dependence of the nucleon self- 
energy as found in [4]. The latter representation is not only more consistent 
with the approximation scheme of the current DBHF approach but also works 
correctly at the Hartree-Fock level using the pseudo- vector pion exchange. This 
minimal requirement, namely that the complete procedure is able to reproduce 
the correct Hartree-Fock results, was never verified for other representations 
and indeed, the 'conventional' pv representation fails in this respect. Hence, a 
strong momentum dependence observed for the nucleon self-energy [4] appears 
to be the artifact of a misrepresentation of the pion exchange potential in the 
previously used 'conventional' pv representations [2,4,15]. 

In this paper we want to go beyond the investigation of Ref. [5] . There we tried 
to determine the range of uncertainty inherent in the present Brueckner ap- 
proach by the consideration of the two limiting cases, namely the pseudo-scalar 
and the complete pseudo- vector representation of the T-matrix. One has, how- 
ever, some additional 'leading order' information on the Lorentz structure of 
the T-matrix, which is given by the Born part, i.e. the bare nucleon-nucleon 
potential. If we want to reproduce the analytically known Hartree-Fock re- 
sults [16] for the complete set of the six non-strange mesons used in the Bonn 
potential, only a mixed 'ps + pv ' representation can be successful where the 
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different parts of the bare potential are represented separately either as ps or 
pv. It is clear that such a mixed representation is not feasible at the DBHF 
level since we can not disentangle the different meson contributions from the 
full in-medium interaction. To proceed, however, we suggest a subtraction 
method to the in-medium T-matrix which means to subtract the bare interac- 
tion, i.e. the leading order, from the full T-matrix and to treat the bare part - 
the Hartree-Fock level - in the mixed 'ps + pv' representation. The remaining 
ladder kernel or 'subtracted T-matrix', as we will call it in the following, i.e. 
the sum of all higher order exchange graphs, is then represented in different 
ways, either as pure ps or 'complete' pv. Applying these two representation 
schemes we are now able to reduce the range of uncertainty concerning the 
T-matrix representation to a minimum. The remaining uncertainty is inherent 
in the current DBHF approach and can not be removed by standard methods. 
The deviations in the final results are, however, small and - similar to the 
treatment of Ref. [5] - the momentum dependence of the self-energy is found 
to be rather weak. Furthermore we obtain new 'Coester lines' for Bonn A, 
B, C with, compared to previous works [3,4], improved saturation properties. 
Most remarkably is a strong softening of the nuclear equation-of-states. 

The paper is now organized as follows: In section 2 we briefly review the Dirac- 
Brueckner Hartree-Fock approach and present some details on the calculation 
of the in-medium on-shell T-matrix elements. In section 3 we introduce the 
projection technique and discuss the different covariant representations used 
for the on-shell T-matrix. The subtraction method is developed at the end of 
this section. The nucleon self-energy in the medium and the nuclear matter 
bulk properties are then presented elaborately in section 4. At the end we 
summarize and conclude our work. 



2 The relativistic Brueckner approach 

2.1 The coupled set of equations 

In the relativistic Brueckner approach the nucleon inside the nuclear medium 
is viewed as a dressed particle in consequence of its two-body interaction with 
the surrounding nucleons. The in-medium interaction of the nucleons is treated 
in the ladder approximation of the relativistic Bethe-Salpeter equation 



where T denotes the T-matrix. V is the bare nucleon-nucleon interaction. The 
intermediate off-shell nucleons in the scattering equation are described by a 
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two-particle propagator iGG. As usually done, we replace this propagator by 
an effective propagator [17], here the Thompson propagator [18] which allows 
only positive energy nucleons in the intermediate scattering states. In addition 
this propagator fixes also the off-shell behavior of the nucleons. This reduces 
the four-dimensional Bethe-Salpeter equation to a three-dimensional integral 
equation. The Pauli operator Q in the Thompson equation accounts for the 
influence of the medium by the Pauli-principle and projects the intermediate 
scattering states out of the Fermi sea. 

The Green's function G which describes the propagation of dressed nucleons 
in the medium fulfills the Dyson equation 

G = Go + G SG . (2) 

G denotes the free nucleon propagator while the influence of the surrounding 
nucleons is expressed by the nucleon self-energy E. In Brueckner theory this 
self-energy is determined by summing up the interaction with all the nucleons 
inside the Fermi sea 

Z = -ij {Tr[GT\ - GT) . (3) 

F 

The Hartree-Fock form of the self-energy integral is necessary if we use an 
'unphysical' T-matrix, as done in [1]. However, since we will entirely work with 
anti-symmetrized two-nucleon states, our T-matrix is 'physical' and contains 
implicitly 'direct' and 'exchange' contributions. Hence, the Hartree form 

£ = -i J Tr[GT] (4) 

F 

of the self-energy integral is sufficient in our case. We will come back to this 
point in more detail in section 3. Since the three equations (1), (2) and (4) 
are strongly coupled, one has to solve this set of equations self-consistently. 

The Dirac structure of the self-energy in isospin saturated nuclear matter 
follows from translational and rotational invariance, parity conservation and 
time reversal invariance [19]. In the nuclear matter rest frame the self-energy 
has the simple form 

k F ) = T, s (k, k F ) -7 S o (A;, k F ) +7 • kE v (fc, k F ) , (5) 

with k^ being the nucleon four-momentum. The self-energy components de- 
pend as Lorentz scalar functions on the Lorentz invariants k 2 , k ■ j and j 2 , 
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where denotes the baryon current. In the nuclear matter rest frame this cur- 
rent is identical to j M = pS^Q, with p being the nuclear matter density. Hence, 
the Lorentz invariants can be expressed in terms of k , |k| and k F , where k F 
denotes the Fermi momentum, related to the density via p = 2/c f /(37t 2 ). By 
taking the traces in Dirac space as [1,4] 

Vs = \tr\E] , S o = ^r[7 S] , £ v = ^±tr [7 • k £] (6) 

one can calculate the different Lorentz components of the self-energy from the 
self-energy operator (4). 

The presence of the medium leads to effective masses and effective momenta 
m*(k,k F ) = M + ReJ: s (k,k F ) , k* = k^ + Re^k, k F ) (7) 

of the nucleons. Above the Fermi surface the self-energy is generally complex 
due to possible decay of particle states into hole states within the Fermi sea. 
To simplify the self-consistency scheme we neglect this decay process and work 
in the so called 'quasi-particle approximation'. Since we only deal with the real 
part of the self-energy we omit this in the notation from now on. 

Defining reduced effective quantities [4,5] 

K = k V ( X + S v(fc, fcp)) , fh*(k, fc F ) = m*(k, bp) I (1 + E v (fc, fc F )) , (8) 

the Dirac equation in the nuclear matter rest frame can be rewritten as 

[ 7m F m - m*(k, k F )}u{k, k F ) = (9) 

which resembles a quasi-free Dirac equation for dressed nucleons. In general the 
reduced effective mass is density but also momentum dependent. To simplify 
the calculation, however, we fix the effective mass of the nucleon in the nuclear 
matter rest frame at the reference momentum |k| = kp. In this 'reference 
spectrum approximation' [20] the reduced effective mass m* F = m*(|k| = k F ) 
works as a self-consistency parameter in the current DBHF approach. All 
equations are iterated until rh* F is converged to a fixed value. The 'reference 
spectrum approximation' implies that the self-energy itself is rather weakly 
momentum dependent. At the end of the calculation one has to verify the 
consistency of the assumption ~ S(|k| = k F ) with the outcome of the 

iteration procedure. 
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Utilizing the different approximations discussed above the positive-energy in- 
medium nucleon spinor are given as 



u\{k, k F ) 



\ 



E*(k) +m* F 
2ml 



2A|k| 



XX 



(10) 



\ £*(k)+m* / 



where E*(k) = ^k 2 + m*p. xx above denotes a two-component Pauli spinor 
with A = ±| and the normalization of the Dirac spinor is u\(k, kp)u\(k, kp) = 
1. Since the in- medium nucleon spinor contains the reduced effective mass the 
matrix elements of the bare nucleon- nucleon interaction are density dependent. 
This density effect does not appear in non-relativistic Brueckner calculations. 
It is believed that it is the main reason for the success of the DBHF approach 
in describing the saturation of nuclear matter [21]. 



2.2 The in-medium T-matrix 



We apply the relativistic Thompson equation [18] to solve the scattering prob- 
lem of two nucleons in the nuclear medium. In the two-particle center of mass 
(cm.) frame - the natural frame for studying the two-particle scattering pro- 
cess - this Thompson equation can be written as [2,4] 



T(p,q,a;)| c . m . = y(p,q) (11) 

r d 3 k , w mf Q(k, x) 

+ / 7 — ^V p,kWi = ^ — ' T(k, q, x) , 

J (2vr) 3 VH ' V 2 (k)2£*(q) -2E*(k) + ie 

where q = (qi — q2)/2 is the relative three-momentum of the initial state 
while k and p are the relative three-momenta of the intermediate and final 
states, respectively. The total four-momentum of the two-nucleon system is 
P* — Qi + <?2> which in the cm. frame becomes P* = (P *,0). P * = = 



2E*(q) = 2yq 2 + m*p is the starting energy in (11). If qi and q 2 are nuclear 
matter rest frame momenta of the nucleons in the initial state, the boost- 
velocity u into the cm. frame is given by 

u = P/V§* + P 2 , (12) 



with the total three-momentum and the invariant mass P = qi + q 2 and 
s* = (£*(qi) + E*(q 2 )) 2 - P 2 , respectively. In Eq. (12) x denotes the set of 
additional parameters x = {k F ,m* F , |u|} on which the T-matrix depends. 

Applying standard techniques as explained in detail by Erkelenz [8] we solve 
the Thompson equation in the cm. frame and calculate the plane- wave helicity 
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matrix elements of the T-matrix. To determine the self-energy only positive- 
energy T-matrix elements at on-shell points |p| = |q| are necessary since in 
(4) we use instead of the full nucleon propagator the Dirac propagator [1,19] 

G D (q) = [7m<T + rh* F }2m5(q* 2 - m^)0(g o *)0(A; F - |q|) . (13) 

Here q denotes the nuclear matter rest frame momentum of the nucleon in the 
Fermi sea. This momentum is on-mass shell, therefore only elastic scattering 
amplitudes, i.e. on-shell T-matrix elements, contribute to the nucleon self- 
energy. Due to the G-functions in the propagator only positive energy nucleons 
are allowed in the intermediate scattering state. Hence the subspace of negative 
energy states is omitted in the current Brueckner approach. In this way we 
avoid the delicate problem of infinities in the theory which would generally 
appear if we would include contributions from negative energy nucleons in the 
Dirac sea [1,6,22]. 

In the on-shell case parity and spin conservation demand that only five of 
sixteen possible positive-energy helicity matrix elements are linearly indepen- 
dent [17]. The five matrix elements are determined explicitly via the |JMLS > 
scheme. After a partial- wave projection onto the | JMLS > states, the Thomp- 
son equation reduces to a partially decoupled set of one-dimensional integral 
equations over the relative momentum |k|. To accomplish such a reduction an 
angle- averaged Pauli-operator Q is used instead of the full Pauli-operator Q. 
Since the Fermi sphere is deformed to a Fermi ellipsoid in the two-nucleon cm. 
frame, Q has to be evaluated for such a Fermi ellipsoid. The explicit expression 
for Q can be found in Refs. [2,4,23]. Finally, the integral equations are solved 
numerically by the matrix inversion technique of Haftel and Tabakin [24] . 

Since the two-nucleon states are two-fermion states, we actually have to eval- 
uate the fully anti-symmetrized matrix elements. Only these matrix elements 
are physically meaningful. Anti-symmetrization is achieved by restoring the 
total isospin of the two-nucleon system (I = 0, 1) with the help of the standard 
selection rule 

(_l)L+s+i = _i _ ( 14 ) 

The five on-shell plane-wave helicity matrix elements for definite isospin are 
finally calculated from the five partial-wave helicity matrix elements obtained 
in the |JMLS > scheme by summing over the total angular momentum J as 

< P A' 1 A;|T I (x)|qA 1 A 2 >= ]T (^p) 4aW < pA' 1 A;|T J ' I (x)|qA 1 A 2 > . 

(15) 
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9 is the scattering angle between q and p, with |p| = |q|, while A = Ai — A2 
and A = A 1 — A 2 . The reduced rotation matrices d 3 x , x (0) are those of Rose [25]. 



3 Covariant representations and the nucleon self-energy 

3. 1 Pseudo-scalar representation 

To use the trace formulas, Eqs. (6), for the calculation of the nucleon self- 
energy components one has to represent the T-matrix in the nuclear matter 
rest frame. Since we determine the T-matrix elements in the two-particle cm. 
frame a representation with covariant operators and Lorentz invariant am- 
plitudes in Dirac space is the most convenient way to Lorentz-transform the 
T-matrix from one frame into another [1]. A set of five linearly independent 
covariants is sufficient for such a T-matrix representation because on-shell only 
five helicity matrix elements appear as solution of the Thompson equation. A 
linearly independent although not unique set of five covariants is given by the 
Fermi covariants 

S = 1 <g> 1, V = Y <g> 7m , T = <g> <v, A = 757^ <g> 757^ P = 75 <g> 75. (16) 

Using this special set - the so called 'pseudo-scalar choice' - the on-shell T- 
matrix for definite isospin I can be represented covariantly as [1] 

T I (\p\,9,x)^F I s (\p\,0,x)S + F 1 ,(\p\,0,x)V + F 1 r (\p\,0,x)T 

+ F I A (\p\,e,x)A + Fl(\p\,9,x)P . (17) 

Here p and 9 denote the relative three-momentum and the scattering angle 
between the scattered nucleons in the cm. frame, respectively. In addition, 
the five Lorentz invariant amplitudes Ff(\p\,0,x) with % = {S,V, T,A,P} 
depend also on x = {kp, m* F , |u|}. We want to stress that the dependence 
of the amplitudes on |p|, 9 and x can be re-expressed in terms of Lorentz 
scalar quantities as explained in detail in [1]. Hence, the representation (17) is 
indeed fully covariant, only the calculation of the Lorentz invariant amplitudes 
is most easily performed in the two-particle cm. frame. 

We determine the Lorentz invariant amplitudes Fj(\p\,9,x) by taking plane- 
wave helicity matrix elements of (17) in the cm. frame. Using as abbreviation 
for the Fermi covariants the operators Tj, with Tj = {S,V, T, A, P}, the on- 
shell matrix elements (|p| = |q|), i.e. the solution of Eq. (15), and the invariant 
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amplitudes are related by 

< P A' 1 A;|T I (x)|qA 1 A 2 >=^<pA' 1 A;|r,|qA 1 A 2 > J F; I (|p|^, a :) . (18) 

i 



Since the Fermi covariants are linearly independent, equation (18) can be 
inverted to determine the unknown amplitudes Fj. The details of this inver- 
sion, especially the treatment of the kinematical singularities at 9 = and 
9 = 7r, which appear likewise in the matrix elements d^, x (9), Eq. (15) and 
< pA' 1 A 2 |rj|qAiA 2 >, are explained in [1,26]. 

Using physical plane-wave helicity matrix elements in (18) the Lorentz invari- 
ant amplitudes F/(|p], 9, x) fulfill a specific anti-symmetry relation which is 
given by the Fierz transformation 

Fl(\p\,n-8,x) = (-l) I Fj(\p\,8,x)r ji , (19) 



where the Fierz matrix T is the matrix given in Eq. (28). Due to this anti- 
symmetry relation the 'direct' representation (17) for the T- matrix is sufficient 
and one can calculate the nucleon self-energy via the Hartree integral (4). 
Using explicitly momenta and Dirac indices, the self-energy reads 

Z af s(k,k F ) = -if Al_G^(q)T(\p\Ax) a ^ T . (20) 



Here k and q denote the incoming and outgoing momenta of the two elastically 
scattered nucleons in the nuclear matter rest frame. The total energy of the 
two nucleons is s* = (E*(k) + E*(q)) 2 — P 2 and the total three-momentum, 
which defines the boost (12) into the two-particle cm. frame, is P = k+q. The 
relative moment um in the c m. frame for which we determine the T-matrix is 
given by |p| = \Js* /A — fn*p. Since we only calculate the Hartree integral, the 
scattering angle is fixed to 9 = 0. 



Applying the covariant representation (17) for the on-shell T-matrix the nu- 
cleon self-energy in isospin saturated nuclear matter is evaluated to be [4] 



d 3 q 9(k F 



(2tt) 3 E*(q) 



(21) 



where the isospin averaged amplitudes are defined as 

Fiflpl^x) := \ [^ I=0 (|p|,0,x) + 3^ I=1 (|p|,0,a;) 



(22) 
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In the self-energy integral (21) only the anti-symmetrized scalar and vector 
amplitudes F$ and Fy for scattering angle 6 = (Hartree) contribute. This is 
not true if we use 'unphysical' helicity matrix elements as done by Horowitz 
and Serot [1]. If we neglect the selection rule (14), we have to determine 
explicitly the Hartree and the Fock contribution to the self-energy via 

r d 4 q 

Hafi(k, k F ) = -i / —— z G° a (q)[T(\p\,0,x) aa ,p T -T(\p\,ir,x) aajT p] .(23) 



An equivalent procedure is to use the Hartree form (20) but with a T-matrix 
representation which explicitly contains 'direct' and 'exchange' terms. This is 
done by calculating 'unphysical' Lorentz invariant amplitudes Fj for scatter- 
ing angles and n — 6, using Eq. (18) with non-antisymmetrized plane-wave 
helicity matrix elements. Then one defines interchanged Fermi covariants as 
[26] 

S = SS , V = SV , T = ST , A = SA , P = SP , (24) 



where the operator S, being the interchange covariant of S, exchanges the 
Dirac indices of particle 1 and 2, i.e. Su(l) a u(2) T = u(l) T u(2) cr . The on-shell 
T-matrix is finally decomposed to 

T l (\pie,x) = T 1 ' D (\p\,e,x)-T 1 ' x (\ P \,e,x) , (25) 



where the 'direct' part of the T-matrix is defined as 



T I ' D (|p|,^x)=[F s I (|p|,^a:)S + ^(|p|^,x)V + ^(|p|^,a:)T 
+ F 1 A (\p\,0,x)A + Fl(\p\,6,x)P , 



(26) 



while the 'exchange' part is given as 



T l ^(\p\,e,x) = (-l) l+1 [F^\p\,n-e,x)S + F^\p\,n-e,x)V 

+F^(\p\,ir-9,x)f + Fl(\p\,ir-e,x)A + Fl(\p\,ir-9,x)P} . (27) 



Since the interchanged and original Fermi covariants are related by the Fierz 
transformation T [26] via 
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the 'exchange' contribution of the T-matrix is identical to 

r^flpl, 0, x) = (-i) I+1 £ ^(| P |, vr - 6i, x^r, 



(28) 



(29) 



Above only the Fermi covariants r, appear and one can rewrite the T-matrix 
representation (25) in terms of a 'direct' representation as in (17). The fully 
anti-symmetrized Lorentz invariant amplitudes are therefore related to the 
non-antisymmetrized amplitudes via the identity 



Fl((\p\,e,x)=Fl((\p\,9,x) 



(30) 



These anti-symmetrized amplitudes respect the relation (19) and are identical 
to the amplitudes which we obtain when we use 'physical' plane-wave helicity 
matrix elements in (18) from the beginning. Nevertheless when employing the 
T-matrix representation (25) with the non-antisymmetrized amplitudes one 
obtains for the self-energy the expression given in Ref. [1] 

*p) = / [&* K - F * + 2 ^v + + Ffj 

(31) 



+rh* F l afS (4F S D - F s x - AF* - 12F* + 4F* - F, 



where the isospin-averaged non-antisymmetrized amplitudes are defined as 
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~ 2 


;F/= (|p|,0,x) + 3F/= 1 ( 


\p\,0,x) 




F* 
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~ 2 


-Ft°(\p\,ir,x) + 3if 


=1 (|p|,7r,a;) 



(32) 



Due to relation (30) the nucleon self-energy defined via the integrals (21) or 
(31) is, of course, identical. 

In Fig. 1 we show the result of a self-consistent DBHF calculation for the 
nucleon self-energy in nuclear matter applying as representation for the on- 
shell T-matrix the ps representation (17). As bare interaction we have used 
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the Bonn A potential [12] and, for comparison, the a-u model potential [29] 
which was originally used by Horowitz and Serot [1]. The density is chosen to 
be the empirical saturation density of nuclear matter with a Fermi momentum 
of k F = 1.34/ra -1 . As already discussed in Ref. [4], we see a pronounced mo- 
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Fig. 1. Momentum dependence of the DBHF nucleon self-energy components in 
nuclear matter at k F = l.Mfm^ 1 using as bare nucleon-nucleon potential Bonn A 
(solid) and the a-iv model potential (dashed). For the T-matrix the ps representation 
(17) is applied. 

mentum dependence of the nucleon self-energy components with the full Bonn 
A while in the case of the a-ou model potential the dependence on the mo- 
mentum is rather weak. A strong momentum dependence questions, of course, 
the validity of the 'reference spectrum approximation' used in the present self- 
consistency scheme. Furthermore, such a strong momentum dependence leads 
to unphysical results deep inside the Fermi sea since the effective mass drops 
to values which are close to zero. Therefore in Ref. [5] the strong momentum 
dependence of the self-energy was studied in more detail and found to orig- 
inate mainly from the one-pion exchange contribution to the self-energy. To 
illustrate this aspect, in Fig. 2 the result of a non-selfconsistent Hartree-Fock 
(HF) calculation are shown. The HF nucleon self-energy is defined via the 
integral 



Z HF = -iJ Tr[G D V] 



(33) 



As in the case of the full DBHF calculation we determine at first the matrix 
elements of V and apply afterwards the ps representation (17) for the matrix 
elements. For better comparison to Fig. 1 we have fixed the Fermi momentum 
again at k F = 1.34/m _1 while for the effective mass, necessary to calculate 
the propagator G D and the dressed potential V, we have taken a fixed value of 
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m* F = 500MeV. With the full Bonn potential the Hartree-Fock nucleon self- 



1000 




k [MeV] 

Fig. 2. Momentum dependence of the Hartree-Fock nucleon self-energy components 
in nuclear matter at kp = 1.34/m -1 . The reduced effective mass of the nucleon is 
fixed at m* F = 500MeV. As nucleon-nucleon interaction the Bonn A potential with 
(solid) and without (dashed) pion exchange part is used. For the bare interaction 
the ps representation (17) is applied. 



energy is similar to the result of the self-consistent DBHF calculation. On the 
other hand applying Bonn A without the pion exchange part, the Hartree-Fock 
nucleon self-energy is weakly momentum dependent, as previously observed in 
the DBHF calculation with the a-ou model potential. This shows that the pion 
exchange part within the Bonn potential is responsible for the strong mo- 
mentum dependence in the nucleon self-energy. That the non-selfconsistent 
Hartree-Fock nucleon self-energy is also similar to the DBHF nucleon self- 
energy in the case of the Bonn A potential indicates that the summation of 
the higher order ladder diagrams of the pion exchange has almost no effect 
on the form of the nucleon self-energy. It appears that the one-pion exchange 
dominates the whole momentum dependence in the nucleon self-energy. Al- 
though the HF and the DBHF self-energies are similar on the scale of some 
100 MeV, one should, however, not conclude that the higher order correlations 
in the T-matrix are of minor importance. Since physical observables, like the 
single particle potential or the equation-of-state result from the cancelation 
of the large scalar and vector fields, they react extremely sensitive on small 
fluctuations on the scale of the self-energy. This is reflected, e.g., by the values 
of the binding energies per particle which is (at fcp = 1.34/m" 1 for full Bonn 
A) E/A = -16.9MeV for the DBHF and E/A = +19.3MeV for the pure HF 
calculation. We will now consider the role of the pion in more detail. 
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3.2 Conventional pseudo-vector representation 



The discussion of the pion-exchange contribution to the nucleon self-energy 
is intimately connected to an ambiguity of the T-matrix representation, as 
it was pointed out in [26]. The set of five covariants used to represent the 
on-shell T-matrix is not uniquely defined when one works exclusively in the 
subspace of positive energy states [22]. Obviously, various alternative sets of 
five linearly independent covariants exist which all can reproduce, like the 
Fermi covariants, the five on-shell helicity matrix elements of the T-matrix. 
For example, the pseudo-vector covariant 

PV = ^£#75 ® 4^75 (34) 



with q* and p* being the initial and final momenta of the nucleons, leads 
to identical on-shell helicity matrix elements as the pseudo-scalar covariant 
P = 75 <E> 75. On-shell the nucleon spinors fulfill the quasi- free Dirac equation 
(9), with fh* F fixed at a reference point, and therefore it holds 

ft*- fa* 

U (P> M - ; l5u{q, k F ) = u{p, fc F )75M(g, fc F ) • (35) 



Thus, if we replace the pseudo-scalar covariant P in the T-matrix representa- 
tion (17) by the pseudo- vector covariant PV and perform the inversion of Eq. 
(18), the calculated Lorentz invariant amplitude Fpy will be identical to the 
previously calculated amplitude F P . Hence, the representation of the on-shell 
T-matrix is ambiguous in the detailed form of the covariant operators one can 
use which is crucial for the description of the one-pion exchange [5] . The 7rNN 
vertex in the OBE potentials, e.g. Bonn, is usually treated by a pseudo-vector 
coupling. In [27] in was shown that the 7rNN coupling is by less than 5% of 
pseudo-scalar nature at the on-shell point. There are several arguments which 
in addition strongly support a pseudo- vector vertex. First of all a PV vertex is 
consistent with soft pion theorems based on chiral symmetry considerations of 
QCD [14]. Secondly, the PV vertex suppresses the coupling to negative energy 
states due to the on-shell relation 

ft* — ft* 

v(p,k F ) * 1 "f 5 u{q,k F ) = , (36) 



with v(p, k F ) a negative energy spinor. In [26] it was, e.g., shown that the one- 
pion exchange contribution to the nuclear optical potential tends to increase 
drastically at low momenta if the 7rNN vertex is treated as pseudo-scalar. 
One reason for this behavior is the strong coupling to negative energy states 
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which is not apparent in non-relativistic approaches. A pseudo-vector vertex 
is more consistent with the approximation scheme of the conventional Brueck- 
ner scheme where one neglects the negative energy states completely. It also 
strongly suppresses the pion contribution in particular at low energies which is 
more in accordance with the empirical knowledge from proton-nucleus scatter- 
ing. One should expect therefore that in the DBHF approach a pseudo-vector 
coupling also drastically reduces the influence of the pion on the nucleon self- 
energy in the medium [1,2]. 

To account for the pseudo- vector nature of the pion exchange in the T-matrix, 
in the past the pseudo-scalar covariant P was simply replaced by the pseudo- 
vector covariant PV [2,4,28,10]. If one uses, however, the 'direct' ps represen- 
tation of the T-matrix (17) and changes to a pv representation, the pseudo- 
vector amplitude F PV = Fp will not contribute to the self-energy (21) since it 
does not appear in the Hartree integral. To suppress the pion contribution one 
needs to find a different pv representation of the T-matrix where the pseudo- 
vector covariant occurs in the Fock or 'exchange' part of the self-energy. One 
possible pv representation is given by the ps representation (25) with P and 
P replaced by PV and PV, respectively. The interchanged covariant PV is 
thereby defined as in [26] by applying the operator S, (24), to the covariant 
PV, (34), and by interchanging the final momenta p\ and p\ in Eq. (34). Per- 
forming this replacement the pseudo-vector amplitude Fp\ = Fp contributes 
to the nucleon self-energy integral with a weight factor of [5] 

TrU* + rh* F )PV] = -( / k* + m* F )(^ F --) . (37) 



For comparison, with the ps representation (25) the weight factor of the 
pseudo-scalar amplitude Fp in the self-energy integral (31) is 

Tr[(jr + rh* F )P] = -(jr-rh* F ) . (38) 



Using the pv representation of the T-matrix as discussed above the nucleon 
self-energy becomes [2,4] 



/d < 



d 3 q Q{k F - |q| 



(2tt) 3 4£*(q) 



( fca/3 /tat 



2q*^k^-q*^) 
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X 
'PV 



+A 



<X0 



AFy - F# + 2F* + 2F£ - 



x 



4mf 



PV 



(39) 
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However, this self-energy appears to be unphysical since the amplitudes F { ' 
used in the integral are calculated from non-antisymmetrized helicity matrix 
elements of the T-matrix. It is by no means clear if the unphysical contribu- 
tions do cancel as they do in the case when we use the ps representation of the 
T-matrix, i.e. see integrals (31) and (21). That the current pv representation 
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Fig. 3. Momentum dependence of the DBHF nucleon self-energy components in 
nuclear matter at A;f = 1.34/m _1 using the Bonn A potential. For the T-matrix the 
ps representation (Eq. (17), solid) and the 'unphysical' pv representation (as used 
in (39), dashed) are applied. 



of the T-matrix is indeed unphysical can be seen in Fig. 3 where the nucleon 
self-energy components, evaluated from the integral (39), are shown. The self- 
energy components continue to be strongly momentum dependent. Further- 
more, the vector component T, v shows an unphysical asymptotic behavior. 
Instead of dropping to zero its value increases with increasing momentum of 
the nucleon. This demonstrates that the pv representation of the T-matrix 
with non-antisymmetrized Lorentz invariant amplitudes Fj, as used in Eq. 
(39), is not useful. 

To circumvent the problem of unphysical contribution to the nucleon self- 
energy one should start from a different ps representation of the T-matrix 
where one uses anti-symmetrized amplitudes Fj(\p\,9,x) and F/(]p|, ir — 9, x) 
instead of unphysical amplitudes -F/(|p|, 9, x) and F/(|p|, ir— 9, x), respectively. 
One possible ps representation of the T-matrix is, for example, given by [4] 

T\\p\,9,x)=T 1 > D (\pl9,x)-T 1 > x (\pl9,x) . (40) 



where the 'direct' and 'exchange' parts of the T-matrix are defined as 
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T 1 ' D (\p\,9,x)^[Fl(\p\,9,x)S + F^(\p\,9,x)Y + F^(\p\,9,x)T 



+ F{(\p\,e,x)A + Fl(\p\,e,x)P 



(41) 



and 



T^ x (\p\,9,x)^(-iy +1 -[Fi(\p\,7r-0,x)S + F^(\p\,7r-9,x)Y 

+F*(\p\, 7T - 9, X)f + F{(\p\, 71 - 0, x)A + i^(j P |,7T - 0, x)P] . (42) 

Due to the anti-symmetry relation (19) for the Lorentz invariant amplitudes 
Fl(\p\,9,x) we have the identity 



T^(\ P \,e,x) = -T^(\p\,e,x) 



(43) 



and with the normalization factors | in (41) and (42) this leads to identical 
results for the self-energies as in the former cases, i.e. Eqs. (21) and (31). 

If we replace now in (41) and (42) the covariants P, P by PV, PV, respectively, 
we arrive at the 'conventional' pv representation as applied in Refs. [2,4,15]. 
The nucleon self-energy components calculated with this 'conventional' pv 
representation of the T-matrix are shown in Fig. 4. The behavior of the nucleon 
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Fig. 4. Momentum dependence of the DBHF nucleon self-energy components in 
nuclear matter at /cf = 1.34/m -1 using the Bonn A potential. For the T-matrix 
the ps representation (Eq. (17), solid) and the 'conventional' pv representation (see 
text after Eq. (43), dashed) are applied. 

self-energy components is similar to the case with the ps representation for 
the T-matrix, i.e. the momentum dependence is still pronounced, although 
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the absolute magnitude of the fields is reduced. As discussed in Ref. [5] this 
is due to the fact that the pion contribution is yet not correctly treated as 
pseudo- vector. 



3.3 Complete pseudo-vector representation 



To suppress the undesirable pseudo-scalar contribution of the pion to the nu- 
cleon self-energy we have to determine a different pv representation of the 
T-matrix. As starting point we should consider the Hartree-Fock level in more 
detail since the correct pv representation of the T-matrix should reproduce 
the HF nucleon self-energy when the pseudo- vector pion exchange potential is 
employed as bare nucleon-nucleon interaction. The 'conventional' pv represen- 
tation discussed in the last section does not fulfill this minimal requirement. 
The reason is simply that due to the Fierz transformation (28) all Fermi covari- 
ants still contain pseudo-scalar contribution, i.e. the covariants (S,V,T,A,P) 
can be re-expressed in terms of (S, V, T, A, P). Hence, in the 'conventional' pv 
representation not all possible pseudo-scalar contributions are replaced by a 
pseudo- vector contribution when one simply replaces P with PV in Eq. (42). 
To obtain a 'complete' pv representation the identities 



i(T + T) = S + S + P + P (44) 
V + V = S + S-P-P (45) 

are actually very helpful. Since Tjon and Wallace already addressed this point 
in Ref. [26] we will follow now partially their notation. To obtain the 'complete' 
pv representation these authors started from a 'symmetrized' ps representation 
of the form 



T I (|p|,^,a:) = / 1 I (|p|,^, a :)(S-S) + / 2 I (|p|,^, a :)-(T + T) 

-/3 I (|p|,^, a :)(A-A) + / 4 I (| P |,^, a ;)(V+V) + / 5 I (|p|,^, a :)(P-P) (46) 

for the T-matrix. Again, this ps representation is equivalent to the ps repre- 
sentation (17). Due to the Fierz transformation (28) the five amplitudes /* 
are related to the amplitudes Fj by a linear transformation 
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(47) 



The anti-symmetry relation (19) converts to the much simpler phase relation 
[26] 



//(|pU-m) = (-) i+ 7/(IpI,M) 



(48) 



where i runs from 1 to 5. Applying the operator identities (45) the 'sym- 
metrized' ps representation (46) can be rewritten as 



T\\p\,6,x)=gl(\p\,6,x)S-gl(\p\,6,x)S + g l A (\p\,e,x)(A-A) 



+ g I P (\p\,9,x)P-gl(\p\,e,x)P 
where the new amplitudes gj are defined as 



(49) 
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(50) 



The relation between the amplitudes gj and Fj is given by the matrix product 
of the linear transformations (47) and (50), i.e. 
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(51) 



Due to the linear relations between the amplitudes Fj, fj and gj, all three ps 
representations (17), (46) and (49) of the T-matrix lead to identical results for 
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the nucleoli self-energy. Using the representation (49) the self-energy integral 
reads simply 



' **) = J (2tt)3 4 ff*( q ) aP ( 93 ~ 9 * 9A ~ 9p) 

+ ^(-9s + 2 9 a + 9p)} , (52) 

where again g { denote isospin averaged amplitudes, evaluated at relative mo- 
menta |p| and scattering angle 9 = in the cm. frame. 

If we replace in (49) the covariants P, P by the pseudo- vector covariants PV, 
PV, respectively, we arrive at the 'complete' pv representation [26] 



T\\p\,9,x)^gl(\p\,9,x)S-gl(\p\,9,x)S + g 1 A (\p\,9,x)(A-A) 

+ 2 PV (|p|,M)PV-4 v (|p|,MPV , (53) 

with gpy(9) and g l ~{9) being identical to g\>{9) and gp(9), respectively. As 
shown in Ref . [5] , this representation is able to reproduce the Hartree-Fock re- 
sults for the nucleon self-energy when we use the pseudo- vector pion exchange 
potential as bare nucleon-nucleon interaction. The self-energy integral using 
this 'complete' pv representation of the T-matrix is given by 



T, af3 (k, fc F ) = J 



(2tt) 3 AE*(q) \ [/Cal3 MaP ' A ™* 2 9pv 
+rh* F l afS Ag s -g§ + %A 

+ ff*p ~9s + 2 9A- 



Amf 

*2 • 



Amy 
Afnf 



PV 



• (54) 



In Fig. 5 we show the non-selfconsistent Hartree-Fock nucleon self-energy com- 
ponents in nuclear matter using only the pion exchange potential as bare in- 
teraction between the nucleons. In the upper panel the ps representation as 
well as the 'conventional' pv representation are used while the lower panel 
shows the results with the 'complete' pv representation of the interaction. 
With the ps representation and the 'complete' pv representation one can re- 
produce the analytical results for the Hartree-Fock nucleon self-energies if we 
use either the pseudo-scalar or the pseudo- vector 7rNN vertex for the pion ex- 
change potential, respectively. Hence, the 'complete' pv representation is the 
correct pseudo-vector representation for the interaction. On the other hand 
the 'conventional' pv representation leads to wrong results for the HF nucleon 
self-energy. In this representation the pion is not completely treated as pseudo- 
vector. Fig. 5 also shows the influence of the single-pion exchange potential 
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Fig. 5. Momentum dependence of the non-selfconsistent Hartree-Fock nucleon 
self-energy components in nuclear matter at kp = 1.34/m -1 . As bare nu- 
cleon-nucleon interaction the single-pion exchange potential is used. In the upper 
panel the ps representation (Eq. (17, solid) and the 'conventional' pv representation 
(see text after Eq. (43), dashed) of the interaction are used. The lower panel shows 
the result using the 'complete' pv representation (Eq. (53), solid) of the interaction. 

to the nucleon self-energy. Only when we use the 'complete' pv representation 
the contribution of the pion to the nucleon self-energy is weak. In all other 
cases, pseudo-scalar or 'conventional' pseudo-vector, the influence of the pion 
is extremely strong, i.e. the contribution to the self-energy is at least one or 
two orders of magnitude larger than in the complete pseudo- vector case. 



In Fig. 6 we present the full self-consistent DBHF calculation with the 'com- 
pete' pv representation of the T-matrix [5]. The DBHF nucleon self-energy 
components are indeed weakly momentum dependent. The single-pion ex- 
change contribution to the interaction, which was previously dominating at 
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Fig. 6. Momentum dependence of the DBHF nucleon self-energy components in 
nuclear matter at = 1.34/m -1 using Bonn A as bare nucleon-nucleon interaction. 
For the T-matrix the ps representation (Eq. (17), solid) and the 'complete' pv 
representation (Eq. (53), dashed) are applied. 

low nucleon momenta, is now strongly suppressed. Consequently the result 
within the 'complete' pv representation using the Bonn A potential resembles 
the result within the a-ui model potential, see Fig. 1 where the ps representa- 
tion was used. To suppress the pion contribution to the in-medium T-matrix 
a correct pseudo- vector like covariant representation is essential for the calcu- 
lation of the nucleon self-energy in nuclear matter. As it is necessary for the 
whole calculation scheme, the weak momentum dependence of the nucleon 
self-energy is also in accordance with the 'reference spectrum approximation' 
used in the calculation. The current DBHF approach therefore appears to be 
self-consistent. 

3.4 Covariant representations of the subtracted T-matrix 

The 'complete' pv representation successfully reproduces the HF nucleon self- 
energy in the case of the pion exchange. Hence, this representation is at the 
moment the 'best' representation of the on-shell T-matrix which is accordance 
with the pseudo-vector nature of the pion exchange potential. However, as al- 
ready pointed out in [5], the 'complete' pv representation fails to reproduce 
the HF nucleon self-energy if other meson exchange potentials are applied 
as bare interaction. This is demonstrated in Fig. 7 where we consider as an 
example the single-omega exchange. While the ps representation of the inter- 
action correctly reproduces the analytical HF nucleon self-energy, the 'com- 
plete' pv representation fails in this respect. In particular the scalar and vector 
self-energies are shifted by about 200 MeV. The failure of the 'complete' pv 
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Fig. 7. Momentum dependence of the non-selfconsistent Hartree-Fock nucleon 
self-energy components in nuclear matter at kp = 1.34 fm~ l using as bare inter- 
action the omega exchange potential. For the potential the ps representation (Eq. 
(17), solid) and the 'complete' pv representation (Eq. (53), dashed) are used. 

representation is understandable since the vector covariant V in the ps repre- 
sentation is partially replaced by a pseudo-vector covariant when we use the 
identity (45) with P and P replaced by PV and PV, respectively. Since the 
'complete' pv representation is not the correct covariant representation of the 
bare interaction, we therefore can not expect that it is correct on the level of 
the full in-medium interaction. 



However, in the last section we have seen that the influence of the pion is 
dominantly given by the single-pion exchange. Hence, it should be reason- 
able to treat the bare interaction and the higher order ladder graphs of the 
meson exchange potential separately. Since the single-meson exchange poten- 
tial is actually known analytically we can represent it covariantly by a mixed 
representation of the form 

V = V pv + V p x . (55) 



Here the n- and 77-meson contributions are treated as pseudo- vector (49) while 
for the (a, cu, p, <5)-meson contributions of the Bonn potential the ps represen- 
tation (46) is applied. The higher order ladder diagrams of the T-matrix 

T Sub = T -V = i I VQGGT = f) f V{iQGGV) n , (56) 

J 71=1 J 

in the following called the subtracted T-matrix, can not be represented cor- 
rectly in a mixed form since we can not disentangle the different meson con- 
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tributions to this part of the full in-medium interaction. The representation 
of the subtracted T-matrix remains therefore ambiguous. However, if the pion 
exchange dominantly contributes to the Hartree-Fock level a ps representa- 
tion of the subtracted T-matrix should be more appropriate because then the 
higher order contributions of other meson exchange potentials are not treated 
incorrectly as pseudo-vector. Thus the most favorable representation of the 
T-matrix is given by the ps representation 

T P = Tl b + V™ + K p ^ s . (57) 
Here the ps representation for Tg ub is determined via the matrix elements 

< pAiAM^aOlqAiAa >:=< pA'^T 1 ^) - \/ I (x)|qA 1 A 2 > , (58) 

with subsequently applying the projection scheme as in Eq. (??). An alterna- 
tive representation of the T-matrix is given by the pv representation 

T PV = TZ + V™ + V ^ S , (59) 

where the subtracted T-matrix is represented by the 'complete' pv represen- 
tation (53). This representation is similar to the 'complete' pv representation 
of the full T-matrix, however, with the advantage that now the pseudo-scalar 
contributions in the bare nucleon-nucleon interaction, e.g. the single-omega 
exchange potential, are represented correctly. In the next section we will use 
both representations, (57) and (59), to study the properties of nuclear mat- 
ter in the DBHF approach. In this way we can determine the influence of the 
higher order ladder graphs to the in-medium interaction in a more quantitative 
way. Furthermore, these two representations set the range of the remaining 
ambiguity concerning the representation of the T-matrix, i.e. after separating 
the leading order contributions. 



4 Results for nuclear matter 

In this section we will present results for the properties of nuclear matter using 
the new approach to the DBHF problem outlined in the previous section. As 
bare nucleon-nucleon potential we employ the one-boson exchange potentials 
Bonn A, B and C [12]. These potentials are based on the exchange of six 
non-strange bosons (ir,rj, p,u,5,a) with masses below 1 GeV. For the pion 
and the eta meson the derivative pseudo- vector coupling is applied. The three 
parameterizations A,B and C of the Bonn potential differ essentially in the 
7rNN form factor and, as a consequence, in the strength of the nuclear tensor 
force. 
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4-1 The nucleon self-energy 
4-1.1 Momentum dependence 

The momentum dependence of the nucleon self-energy at saturation density 
in isospin saturated nuclear matter is shown in Fig. 8. For both T- matrix 
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Fig. 8. Momentum dependence of the nucleon self-energy components in nuclear 
matter at kp = 1.34/m _1 using Bonn A as bare nucleon-nucleon interaction. For 
the T-matrix the subtraction scheme with the ps representation (Eq. (57), solid) 
and the pv representation (Eq. (59), dashed) are applied. 

representations (57) and (59), the self-energy components are rather weakly 
dependent on the nucleon momentum. In addition, they are also almost identi- 
cal, i.e. the difference between a pseudo-scalar or pseudo- vector representation 
of the higher order ladder diagrams is rather small. This demonstrates again 
that the pion indeed contributes mostly to the Hartree-Fock level. Hence the 
ambiguity of the T-matrix representation has only minor influence on the final 
result for the nucleon self-energy in the medium 

However, at larger densities of the nuclear medium the situation changes. In 
Fig. 9 the momentum dependence of the self-energy components in nuclear 
matter at a Fermi momentum of hp = 1.8/m -1 is shown. Now the two repre- 
sentations of the T-matrix lead to significantly different results. In addition, 
the momentum dependence of the nucleon self-energy components is increasing 
for both representations. A strong momentum dependence of the self-energies 
at higher densities was also observed in Ref . [2] . The question how to represent 
the T-matrix is now much more severe than at lower densities. However, the 
'complete' pv representation of the subtracted T-matrix gives a rather unphys- 
ical asymptotic behavior of the nucleon self-energy. The vector component Sy 
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Fig. 9. Momentum dependence of the nucleon self-energy components in nuclear 
matter at kp = 1.8fm~ 1 using Bonn A as bare nucleon-nucleon interaction. For the 
T- matrix the subtraction scheme with the ps representation (Eq. (57), solid) and 
the pv representation (Eq. (59), dashed) are applied. 

not only changes sign but also increases drastically at large nucleon momenta. 
We believe that the misrepresentation of the higher order ladder diagrams 
of the heavy meson exchange potentials which should be treated as pseudo- 
scalar and not as pseudo-vector is responsible for this behavior. Since the ps 
representation of the ladder kernel, i.e. the subtracted T-matrix, still yields 
reasonable results at higher densities this representation should be preferable. 
Since the momentum dependence increases with increasing density, it should 
be included in the future in the self-consistency scheme when predictions at 
high densities are made [11]. 

In Fig. 10 we show the complete dependence of the scalar S5 self-energy and 
the vector component So on momentum and density. As bare interaction the 
Bonn A potential is used again while for the T-matrix the subtraction scheme 
with the ps representation of the ladder kernel, Eq. (57), is applied. As can 
be seen from Fig. 10 the momentum dependence starts to become pronounced 
at densities around p = 2.0po, where po — 0.166/m~ 3 is the empirical satu- 
ration density of nuclear matter. At fixed density the momentum dependence 
is still most pronounced around the corresponding Fermi momentum. On the 
other hand, keeping the nucleon momentum fixed and varying the Fermi mo- 
mentum, we see that the medium dependence, namely the variation of the 
self-energy with the nuclear matter density, is strongest for low energetic nu- 
cleons. This clearly demonstrates the influence of the Pauli-blo eking effect 
which vanishes with increasing relative momentum of the nucleon interacting 
with the particles inside the Fermi sea. 
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Fig. 10. Dependence of the nucleon self-energy component S s (upper panel) and 
So (lower panel) on the nuclear matter density p and the nucleon momentum k. As 
bare interaction the Bonn A potential is used while for the T-matrix the subtraction 
scheme with the ps representation, Eq. (57), is applied. 

4-1-2 Density dependence 

The detailed density dependence of the nucleon self-energy components in 
nuclear matter is presented in Fig. 11. The momentum k of the nucleon is 
thereby fixed at the Fermi momentum fop. The density dependence of the 
nucleon self-energy is quite similar for both representations (57) and (59) of 
the T-matrix. An important difference is, however, the deviation in the vector 
component Ey. As already seen from Fig. 9, the 'complete' pv representation of 
the ladder kernel leads to an unphysical behavior of the spatial Ey component. 

This has a strong influence on the reduced effective mass fh* F of the nucleon, 
Eq. (8). The reduced effective mass generally drops with increasing density 
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Fig. 11. Density dependence of the nucleon self-energy components in nuclear mat- 
ter using Bonn A as bare nucleon-nucleon interaction. For the T-matrix the subtrac- 
tion scheme with the ps representation (Eq. (57), solid) and the pv representation 
(Eq. (59), dashed) are applied. 

as can be seen from Fig. 12 where the reduced effective mass is shown as 
a function of the Fermi momentum. At saturation density, for Bonn A at 




Fig. 12. Density dependence of the reduced effective mass rh F of the nucleon in 
nuclear matter, using Bonn A as bare nucleon-nucleon interaction. For the T-matrix 
the subtraction scheme with the ps representation (Eq. (57), solid) and the pv 
representation (Eq. (59), dashed) are applied. In addition the result of Horowitz 
and Serot (HS) (Ref. [1], solid with circles) and of Brockmann and Machleidt (BM) 
(Ref. [3], dotted) are shown. 

kp = 1.39/m _1 , the reduced effective mass has a value of rh* F ~ 0.65M for 
both T-matrix representations. This agrees with the findings of other groups 
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[1,3] and also with the empirical value determined from the spin-orbit interac- 
tion in finite nuclei [30]. At larger densities, however, the two representations 
for the T-matrix lead to rather different results. At a Fermi momentum of 
hp = 2.0/m~ 1 , which is three times nuclear matter density, the results differ 
by almost a factor of two. The solid curve with circles shown in Fig. 12 is the 
result of the calculation by Horowitz and Serot [1], who used as bare nucleon- 
nucleon interaction the a-u model potential. Since they did not consider the 
pion in their calculation they used the ps representation for the full in-medium 
T-matrix. Fitting the saturation properties of nuclear matter, their result for 
the effective mass is almost identical to our calculation with the ps represen- 
tation (57) of the subtracted T-matrix. The calculation of Brockmann and 
Machleidt [3] is shown as dotted curve in Fig. 12. These authors used instead 
of the projection technique a fit procedure to the single-particle potential. 
They assumed a momentum independent form of the nucleon self-energy and 
thus they could obtain only approximately the effective mass of the nucleon. 
For lower densities their result is, however, quite similar to our findings. Only 
at higher densities important differences occur. In their calculation the effec- 
tive mass seems to saturate much earlier while in our calculation with the ps 
representation of the subtracted T-matrix the effective mass drops to a smaller 
value at larger densities. 

The above analysis indicates again that the ps representation of the remaining 
ladder kernel of the T-matrix is preferable compared to the 'complete' pv rep- 
resentation. Although the two alternatives proposed here yield similar results 
at densities below and up to saturation, the different representations become 
decisive with increasing density. The 'complete' pv representation of the lad- 
der kernel of the T-matrix leads thereby to a partially unphysical high density 
behavior, i.e. a too strongly dropping mass, whereas the ps representation is 
still reasonable. For completeness the various relevant quantities as a function 
of the nuclear matter density are presented in Tab. 1 using the subtraction 
scheme with the ps representation (57) for the subtracted T-matrix. As bare 
nucleon-nucleon interaction the Bonn A potential is used. 

4-2 The equation- of- state of nuclear matter 

In the relativistic Brueckner theory the energy per particle is defined as the 
kinetic plus half the potential energy 

^M = ^E<«A(k)| 7 -k + M + ls(A;)K(k)>|J||-M . (60) 

In Fig. 13 we show the binding energy per particle E/A as a function of the 
density, calculated with Bonn A, B and C. For the T-matrix the subtraction 
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Table 1 

The Fermi- momentum kp, the nuclear matter density p, the binding energy per 
particle E/A, the reduced effective mass m* F and the components of the nucleon 
self-energy (at |k| = kp) for nuclear matter applying the ps representation (57) for 
the subtracted T-matrix. As bare nucleon-nucleon interaction the Bonn A potential 
is used. 



kp 


p 


E/A 


rh* F 


s s 


-S 




[far 1 ] 


[far 3 ] 


[MeV] 


[MeV] 


[MeV] 


[MeV] 




0.7 


0.023 


-7.74 


868.4 


-72.0 


42.1 


-0.0017 


0.8 


0.035 


-8.77 


858.2 


-85.6 


51.1 


-0.0058 


0.9 


0.049 


-9.98 


837.1 


-108.8 


68.0 


-0.0084 


1.0 


0.068 


-11.62 


809.1 


-139.5 


90.6 


-0.0121 


1.1 


0.090 


-13.45 


771.6 


-179.2 


120.9 


-0.0155 


1.2 


0.117 


-14.80 


729.9 


-224.0 


157.3 


-0.0207 


1.3 


0.148 


-15.74 


682.8 


-275.2 


200.3 


-0.0281 


1.35 


0.166 


-16.03 


657.8 


-303.0 


224.2 


-0.0333 


1.4 


0.185 


-16.15 


632.1 


-331.8 


249.5 


-0.0395 


1.5 


0.228 


-15.28 


577.8 


-392.6 


305.8 


-0.0546 


1.6 


0.277 


-12.37 


520.4 


-457.0 


371.1 


-0.0741 


1.7 


0.332 


-6.19 


462.9 


-522.2 


445.8 


-0.0999 


1.8 


0.394 


4.78 


408.4 


-584.6 


530.2 


-0.1325 


1.9 


0.463 


22.22 


358.9 


-641.4 


626.4 


-0.1712 


2.0 


0.540 


48.19 


315.9 


-690.8 


735.6 


-0.2146 



scheme with the ps representation (57) is applied. With Bonn A one can repro- 
duce the empirical saturation point of nuclear matter, shown as shaded region 
in the figure. The other Bonn potentials give less binding energy although the 
saturation density is always close to the empirically known value. 

The result for the binding energy per particle using the two representations 
(57) and (59) for the T-matrix are very similar as can be seen in Fig. 14. 
At saturation density the binding energy is only 0.5 MeV smaller using the 
pseudo- vector representation of the subtracted T-matrix. Thus, the energy per 
particle is not very sensitive on the explicit representation of the subtracted 
T-matrix. As already noticed, on the level of the self-energies, Fig. 5, 9 and 12, 
differences between the two methods occur at higher densities. In particular 
the equation-of-state in the ps representation for the subtracted T-matrix 
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Fig. 13. Binding energy per particle as a function of nuclear matter density. As 
bare nucleon-nucleon interaction the potentials Bonn A, B and C are used. For the 
T-matrix the subtraction scheme with the ps representation (57) is applied. The 
shaded box denotes the empirical saturation region of nuclear matter. 
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Fig. 14. Binding energy per particle as a function of nuclear matter density. As 
bare nucleon-nucleon interaction the Bonn A potential is used. For the T-matrix 
the subtraction scheme with the ps representation (Eq. (57), solid) and the pv 
representation (Eq. (59), dashed) are applied. 

appears to be more stiff at higher densities than with the corresponding pv 
representation. 



To understand this behavior in more detaif and to compare also with other 
calculations we show in Fig. 15 the binding energy per particle as a function of 
nuclear matter density for different scenarios. First of all, applying the projec- 
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Fig. 15. Binding energy per particle as a function of nuclear matter density. As 
bare nucleon-nucleon interaction the Bonn A potential is used. For the T-matrix 
the subtraction scheme with the ps representation (Eq. (57), solid), the full ps 
representation (Eq. (17), dash-dotted) and the 'conventional' pv representation (see 
text after Eq. (43), dashed) are applied. In addition the result of Brockmann and 
Machleidt (BM) (Ref. [3], dotted) is shown. 



tion method of Horowitz and Serot, we can reproduce almost completely the 
result of Brockmann and Machleidt (BM) , Ref. [3], when we use the 'conven- 
tional' pv representation of the T-matrix. (see text after Eq. (43)). The two 
approaches are in detail compared in Tab. 3. Except of the incompressibility, 
which is generally smaller in our calculation, the bulk properties are indeed 
very similar. This result is somewhat surprising since the approach of Ref. [3] 
is completely different from the present one. As already mentioned, in Ref. 
[3] no projection scheme to the T-matrix has been applied but constant, i.e. 
momentum independent, self-energy components have been determined by a 
fit to the single particle potential. On the other hand, the 'conventional' pv 
representation of the T-matrix leads to a pronounced momentum dependence 
of the self-energy components, see Fig. 4. That these two calculations give 
nevertheless similar results for the nuclear matter bulk properties is due to 
the fact that they lead to more or less identical values for the fields and the 
effective mass at the Fermi momentum. This agreement appears, however, to 
be somewhat fortuitous. 



However, as discussed in the previous sections the 'conventional' pv representa- 
tion does not correctly reproduce the contribution of the single-pion exchange 
potential to the nucleon self-energy. On the level of the binding energy one 
can estimate the contribution of the single-pion exchange potential compar- 
ing the result of ps representation for the subtracted T-matrix and the pure 
ps representation of the T-matrix. In the latter approach the nucleons are 
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less bound at small densities but the situation changes around saturation. A 
correct pseudo-vector representation of the pion, as used in the subtraction 
scheme, suppresses this effect. Thus at smaller densities we obtain a larger 
binding, while around saturation density the binding energy is smaller. Com- 
pared to the 'conventional' pv representation or the result of Brockmann and 
Machleidt, the ps representation of the subtracted T-matrix plus a correct 
treatment of the bare interaction leads altogether to more binding at smaller 
and medium densities. The correct treatment of the single-pion exchange po- 
tential is therefore essential in the low density regime of the equation-of-state. 
Except of the pure ps treatment, which is certainly not correct for a realis- 
tic potential like the Bonn potential [10,2,5], the various calculations coincide 
at high densities. This reflects a decreasing relative importance of the pion- 
exchange at high densities. 

-10 i ■ . ■ . ■ 1 



-12 



> -14 

S3 - 16 



-18 



-20 



• A 

O 



B 



I A 



• ps for l Sub 

o 'complete' pv for T Sub 

a BM 



1.2 



1.3 



k F [fm 1 ] 



1.4 



1.5 



Fig. 16. Saturation points of nuclear matter. As bare nucleon-nucleon interaction 
the Bonn potentials A,B and C are used. For the T-matrix the subtraction scheme 
with the ps representation (Eq. (57), filled circles) and the pv representation (Eq. 
(59), open circles) are applied. As open triangles the results of the calculation of 
Brockmann and Machleidt (BM), Ref. [3], are shown. 

The present results are summarized in Fig. 16 where the corresponding satura- 
tion points for the three different versions of the Bonn potential are shown. We 
compare the results with the two representation of the subtracted T-matrix 
with the results of the calculation of Brockmann and Machleidt (BM), Ref. 
[3]. With the improved representation schemes (57) and (59) for the T-matrix 
one obtains new 'Coester-lines' which are left of the original one, i.e. shifted 
towards the empirical region. The refined treatment of the T-matrix repre- 
sentation leads to an enhancement of the binding energy connected with a 
reduced saturation density. As in the previous calculations, Bonn A is still the 
only one which meets the empirical region. However, due to an increased bind- 
ing energy Bonn B is now much closer to empirical region. This observation 
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is consistent with the present treatment. The different types of Bonn poten- 
tials essentially vary in the strength of the nuclear tensor force determined 
by the 7rNN form factor. Bonn A which has the smallest tensor force yields 
the smallest D-state probability of the deuteron and only a pure description 
of the 3 Di phase shift [3,12]. Thus it appears that a refined treatment of the 
pion exchange leads to improved nuclear matter results for the more realistic 
Bonn B potential. Bonn C, however, is still far off the empirical region. 

Furthermore, it can be seen from Fig. 16 and Tab. 2 that the final nuclear 
matter bulk properties depend only moderate on the representation of the 
subtracted T-matrix. In Ref. [5] we tried already to determine the range of 
inherent uncertainty in the relativistic Brueckner approach which is due to 
the ambiguities concerning the representation of the T-matrix discussed in 
Section 3. By the separate treatment of the Born contribution to the T-matrix 
we end up now with a much narrower uncertainty band which is given by the 
ps or complete pv representation of the ladder kernel, i.e. the subtracted T- 
matrix. Over the different types of Bonn interactions the two methods lead 
to a variation of 0.5 MeV in the binding energy, 0.1-0.2 fm~ l in the Fermi 
momentum, and to about 30 MeV concerning the value of the effective mass 
at saturation density. The values for incompressibility are also close in the 
two approaches, i.e. they differ by less than 10 MeV. However, compared to 
the conventional pv representation and especially to Ref. [3], see Tab. 3, the 
kompression moduli are significantly reduced for all three types of interactions 
when we use the new approach to the T-matrix representation, see Tab. 2. 
Within the ps representation of the subtracted T-matrix, Bonn B and C now 
yield very small kompression moduli around K = 150MeV and K = 115MeV, 
respectively. For Bonn A a value of K = 230MeV is obtained. This value 
agrees with the empirical value of the kompression modulus of K — 210 ± 
30MeV [31]. Here Brockmann and Machleidt found much larger values for all 
three Bonn potentials. 
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Table 2 

The Fermi-momentum kp, the binding energy per particle E/A, the reduced ef- 
fective mass rhp and the kompression modulus K at saturation for nuclear matter 
using as bare nucleon-nucleon interaction the Bonn potentials A, B and C. For 
the T-matrix the subtraction scheme with the ps representation (57) and the pv 
representation (59) are applied. 



ps for Tsub 'completely for Ts u b 



k F E/A m* F K k F E/A m* F K 

[fin" 1 ] [MeV] [MeV] [MeV] [fin" 1 ] [MeV] [MeV] [MeV] 



A 


1.39 


-16.15 


637.0 


230 


1.39 


-15.72 


596.0 


220 


B 


1.34 


-14.59 


667.0 


150 


1.33 


-13.99 


634.0 


140 


C 


1.29 


-13.69 


694.0 


110 


1.27 


-13.00 


669.0 


100 



Table 3 

The Fermi-momentum k F , the binding energy per particle E/A, the reduced ef- 
fective mass rhp and the kompression modulus K at saturation for nuclear matter 
using as bare nucleon-nucleon interaction the Bonn potentials A, B and C. For the 
T-matrix the 'conventional' pv representation, see text after Eq. (43), is applied. In 
addition the results of Brockmann and Machleidt (BM), Ref. [3], are presented. 







'conventional' pv 






BM 






k F 


E/A 


— * 

m F 


K 


k F 


E/A 




K 




[fm- 1 ] 


[MeV] 


[MeV] 


[MeV] 


[fm" 1 ] 


[MeV] 


[MeV] 


[MeV] 


A 


1.41 


-15.81 


538.0 


275 


1.40 


-15.59 


564.0 


290 


B 


1.35 


-13.70 


565.0 


195 


1.37 


-13.60 


573.0 


249 


C 


1.30 


-12.31 


585.0 


155 


1.32 


-12.26 


590.0 


185 
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5 Summary 



We have investigated the nuclear matter properties in the relativistic Brueck- 
ner approach. The required representation of the T-matrix by Lorentz in- 
variant amplitudes suffers thereby from on-shell ambiguities concerning the 
pseudo-scalar or pseudo-vector nature of the interaction. We minimized this 
ambiguity by separating the leading order, i.e. the single-meson exchange, from 
the full T-matrix. Actually we represented the contributions stemming from 
the single-meson exchange by taking the pseudo-scalar and the pseudo- vector 
nature of the interaction into account. Up to now, this approach is the only 
one which reproduces the correct results for the T-matrix on the Hartree-Fock 
level. The remaining higher order correlations, i.e. the ladder kernel, are then 
represented either completely as pseudo-scalar or as pseudo-vector. 

This method takes at best the pion contribution to the nucleon self-energy 
into account which, on the other hand, essentially determines the momentum 
dependence of the self-energy. Treating the one-pion exchange with a pseudo- 
vector coupling the momentum dependence is strongly suppressed as it is 
desired from meson phenomenology. This also favors the 'reference spectrum 
approximation' which was a first guess to the momentum dependence of the 
effective mass of the nucleon. Now the results of the calculation are much 
more consistent with the initial assumption of a weakly momentum dependent 
effective mass. 

Furthermore, we have investigated the density dependence of the nucleon self- 
energy. The results for both scenarios for the T-matrix representation are quite 
similar up to saturation density. This holds also for the binding energy. Both 
representations lead to similar saturation properties of nuclear matter which 
indicates that the remaining ambiguity in the representation of the ladder 
kernel, i.e. the higher order correlations, is not too severe at moderate densi- 
ties. However, the two schemes start to differ significantly at higher densities. 
There, the complete pv representation of the ladder kernel, i.e. the subtracted 
T-matrix, leads to an unphysical behavior, whereas the corresponding ps rep- 
resentation still leads to reasonable results. Thus it appears that the higher 
density correlations are best represented adopting the pseudo-scalar represen- 
tation. 

As a major result of our investigation we obtain new 'Coester lines' for the 
various Bonn potentials. Compared to previous treatments these are shifted 
towards the empirically know saturation point. Bonn A is still the only poten- 
tial which really meets the empirical region of saturation, but, with improved 
saturation properties compared to previous treatments. The refined treatment 
of the pion exchange leads on the other hand also to improved results for the 
- from the view of the phase shift analysis - more realistic Bonn B potential. 
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Furthermore, we found that the equation-of-state is strongly softened com- 
pared to previous calculations. Actually with Bonn A we obtain a kompres- 
sion modulus of K ~ 230MeV which is in good agreement with the empirical 
value. 

To summarize our results, we obtained new results for the nuclear matter 
properties within the projection technique employing an new method for the 
T-matrix representation. The final results are at lower densities almost insensi- 
tive on the explicit choice made for the representation. However, at higher den- 
sities, certain differences occur when using different representation schemes. 
We want to stress that the ambiguity in the projection technique is still not 
fully resolved yet. We plan to look on off-shell T-matrix elements in the future 
since off-shell matrix elements of the pseudo-scalar and pseudo-vector covari- 
ants differ significantly. We hope that this might bring more insight on what 
is the correct on-shell representation of the T-matrix. 
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